1. We load the CITE-seq human PBMC data. You can follow the Seurat tutorial or download the analyzed R object.

  1. We run WLR on the RNA data.
ptm <- proc.time()
cbmc_rna <- RunWLR(object = cbmc_rna, n_iters = 1000, wlr_pca_dims = 20, wlr_k = 30, ratio = 0.1, norm_type = "l1")
Run Bootstrapping
Finalize
ptm1 <- proc.time() 
print(x = ptm1 - ptm)
   user  system elapsed 
 42.979  13.618  70.487 
  1. We run WLR on the ADT data.
print(x = ptm1 - ptm)
   user  system elapsed 
 34.798  11.520  58.947 
  1. We fuse the RNA and ADT data. We see that CD4+ and CD8+ T cells are more defined by the ADT data in the surface protein space, while Dendritic cells (DCs) are more defined in the RNA-seq data due to the absence of its surface marker in the ADT data.

  1. We perform UMAP dimensional reduction on the fused data. Additionally, clustering can be performed as well.

LS0tCnRpdGxlOiAiV0xSIFR1dG9yaWFsIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCjEuIFdlIGxvYWQgdGhlIFtDSVRFLXNlcV0oaHR0cHM6Ly93d3cubmF0dXJlLmNvbS9hcnRpY2xlcy9ubWV0aC40MzgwKSBodW1hbiBQQk1DIGRhdGEuIFlvdSBjYW4gZm9sbG93IHRoZSBbU2V1cmF0IHR1dG9yaWFsXShodHRwczovL3NhdGlqYWxhYi5vcmcvc2V1cmF0L3YzLjEvbXVsdGltb2RhbF92aWduZXR0ZS5odG1sKSBvciBkb3dubG9hZCB0aGUgYW5hbHl6ZWQgW1Igb2JqZWN0XShodHRwczovL2RyaXZlLmdvb2dsZS5jb20vZHJpdmUvZm9sZGVycy8xaG12ZjF2SmRHMkRONjBmbDFqekxXbUNNLVNOd1lqalg/dXNwPXNoYXJpbmcpLgoKYGBge3Igc2V0dXAsIGZpZy5oZWlnaHQ9NSwgZmlnLndpZHRoPTEwfQpsaWJyYXJ5KFdMUikKbGlicmFyeShTZXVyYXQpCmNibWMgPC0gcmVhZFJEUyhmaWxlID0gImNibWMuZGF0YS5yZHMiKQpDb21iaW5lUGxvdHMocGxvdHMgPSBsaXN0KERpbVBsb3QoY2JtYywgcmVkdWN0aW9uID0gInVtYXAiLCBncm91cC5ieSA9ICJybmFDbHVzdGVySUQiLCBsYWJlbCA9IFRSVUUpICsgZ2dwbG90Mjo6Z2d0aXRsZSgiUk5BLXNlcSBEYXRhIiksIERpbVBsb3QoY2JtYzEsIHJlZHVjdGlvbiA9ICJ0c25lX2FkdCIsIGdyb3VwLmJ5ID0gInJuYUNsdXN0ZXJJRCIsIGxhYmVsID0gVFJVRSkgKyBnZ3Bsb3QyOjpnZ3RpdGxlKCJBRFQgU3VyZmFjZSBQcm90ZWluIERhdGEiKSkpCmBgYAoKMi4gV2UgcnVuIFdMUiBvbiB0aGUgUk5BIGRhdGEuCmBgYHtyIFJOQX0KZGF0YV9ybmEgPC0gY2JtY0Bhc3NheXMkUk5BQGNvdW50c1tjYm1jQGFzc2F5cyRSTkFAdmFyLmZlYXR1cmVzLF0KY2JtY19ybmEgPC0gSW5pdGlhdGVXTFIocmF3X2RhdGEgPSBkYXRhX3JuYSwgcHJvamVjdF9uYW1lID0gIlJOQSIsIG5vcm1hbGl6YXRpb24gPSAiY29zaW5lIiwgcGNhX2RpbXMgPSAyMCwga2VybmVsID0gImdhdXNzaWFuIiwgayA9IDMwKQpjYm1jX3JuYSA8LSBHZW9tU2tldGNoKG9iamVjdCA9IGNibWNfcm5hLCBnZW9tX3BjYV9kaW1zID0gMjAsIGdlb21fc2l6ZSA9IDEwMDAsIHNrZXRjaF9uX3BjYSA9IDIwLCBza2V0Y2hfayA9IDMwKQpwdG0gPC0gcHJvYy50aW1lKCkKY2JtY19ybmEgPC0gUnVuV0xSKG9iamVjdCA9IGNibWNfcm5hLCBuX2l0ZXJzID0gMTAwMCwgd2xyX3BjYV9kaW1zID0gMjAsIHdscl9rID0gMzAsIHJhdGlvID0gMC4xLCBub3JtX3R5cGUgPSAibDEiKQpwdG0xIDwtIHByb2MudGltZSgpIApwcmludCh4ID0gcHRtMSAtIHB0bSkKYGBgCgozLiBXZSBydW4gV0xSIG9uIHRoZSBBRFQgZGF0YS4KYGBge3IgQURUfQpkYXRhX3BydCA8LSBjYm1jQGFzc2F5cyRBRFRAY291bnRzCmNibWNfcHJ0IDwtIEluaXRpYXRlV0xSKHJhd19kYXRhID0gZGF0YV9wcnQsIHByb2plY3RfbmFtZSA9ICJQcm90ZWluIiwgbm9ybWFsaXphdGlvbiA9ICJjIiwga2VybmVsID0gImdhdXNzaWFuIiwgayA9IDMwLCB2ZXJib3NlID0gRkFMU0UpCmNibWNfcHJ0IDwtIEdlb21Ta2V0Y2gob2JqZWN0ID0gY2JtY19wcnQsIGdlb21fc2l6ZSA9IDEwMDAsIGdlb21fcGNhX2RpbXMgPSAzLCBza2V0Y2hfbl9wY2EgPSAwLCBza2V0Y2hfayA9IDMwKQpwdG0gPC0gcHJvYy50aW1lKCkKY2JtY19wcnQgPC0gUnVuV0xSKG9iamVjdCA9IGNibWNfcHJ0LCBuX2l0ZXJzID0gMTAwMCwgcmF0aW8gPSAwLjEsIHdscl9rID0gMzAsIG5vcm1fdHlwZSA9ICJsMSIpCnB0bTEgPC0gcHJvYy50aW1lKCkgCnByaW50KHggPSBwdG0xIC0gcHRtKQpgYGAKCjQuIFdlIGZ1c2UgdGhlIFJOQSBhbmQgQURUIGRhdGEuIFdlIHNlZSB0aGF0IENENCsgYW5kIENEOCsgVCBjZWxscyBhcmUgbW9yZSBkZWZpbmVkIGJ5IHRoZSBBRFQgZGF0YSBpbiB0aGUgc3VyZmFjZSBwcm90ZWluIHNwYWNlLCB3aGlsZSBEZW5kcml0aWMgY2VsbHMgKERDcykgYXJlIG1vcmUgZGVmaW5lZCBpbiB0aGUgUk5BLXNlcSBkYXRhIGR1ZSB0byB0aGUgYWJzZW5jZSBvZiBpdHMgc3VyZmFjZSBtYXJrZXIgaW4gdGhlIEFEVCBkYXRhLgpgYGB7ciBmdXNlLCBmaWcuaGVpZ2h0PTUsIGZpZy53aWR0aD0xMH0KZnVzZWQgPC0gRnVzZVdMUihjYm1jX3JuYSwgY2JtY19wcnQsIHByb2plY3RfayA9IDEwKQpjYm1jQG1ldGEuZGF0YVsiUk5BX1Njb3JlIl0gPC0gZnVzZWQkZnVzZWRfd2VpZ2h0WyJSTkEiLF0KY2JtY0BtZXRhLmRhdGFbIlByb3RlaW5fU2NvcmUiXSA8LSBmdXNlZCRmdXNlZF93ZWlnaHRbIlByb3RlaW4iLF0KQ29tYmluZVBsb3RzKHBsb3RzID0gbGlzdChWbG5QbG90KG9iamVjdCA9IGNibWMxLCBmZWF0dXJlcyA9ICJSTkFfU2NvcmUiLCBzb3J0ID0gVFJVRSksIFZsblBsb3Qob2JqZWN0ID0gY2JtYzEsIGZlYXR1cmVzID0gIlByb3RlaW5fU2NvcmUiLCBzb3J0ID0gVFJVRSkpKQpgYGAKCjUuIFdlIHBlcmZvcm0gVU1BUCBkaW1lbnNpb25hbCByZWR1Y3Rpb24gb24gdGhlIGZ1c2VkIGRhdGEuIEFkZGl0aW9uYWxseSwgY2x1c3RlcmluZyBjYW4gYmUgcGVyZm9ybWVkIGFzIHdlbGwuCmBgYHtyIHVtYXAsIGZpZy5oZWlnaHQ9NCwgZmlnLndpZHRoPTV9CnVtYXAuZGltcyA8LSB1d290Ojp1bWFwKFggPSBhcy5tYXRyaXgoZnVzZWQkZnVzZWRfZGlzdCksIG5fY29tcG9uZW50cyA9IDIpCmRpbW5hbWVzKHggPSB1bWFwLmRpbXMpIDwtIGRpbW5hbWVzKHggPSBjYm1jQHJlZHVjdGlvbnMkdW1hcEBjZWxsLmVtYmVkZGluZ3MpCmNibWNAcmVkdWN0aW9ucyR1bWFwQGNlbGwuZW1iZWRkaW5ncyA9IHVtYXAuZGltcwpEaW1QbG90KG9iamVjdCA9IGNibWMsIGdyb3VwLmJ5ID0gInJuYUNsdXN0ZXJJRCIsIHJlZHVjdGlvbiA9ICJ1bWFwIiwgbGFiZWwgPSBUUlVFKQpgYGAK